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ABSTRACT 

We calculate the observable signature of a black hole accretion disk with a gap or hole created by a secondary 
black hole embedded in the disk. We find that for an interesting range of parameters of black hole masses 
(~ 10 6 — 10 9 M ), orbital separation (<~ 1 AU to ~ 0.1 pc), and gap width (1-19 disk scale heights), the missing 
thermal emission from a gap manifests itself in an observable decrement in the spectral energy distribution. We 
present observational diagnostics in terms of powerlaw forms that can be fit to line-free regions in AGN spectra 
or in fluxes from sequences of broad filters. Most interestingly, the change in slope in the broken powerlaw is 
almost entirely dependent on the width of gap in the accretion disk, which in turn is uniquely determined by 
mass ratio of the black holes, such that it scales roughly as q 5 / 4 . Thus one can use spectral observations of 
the continuum of bright active galactic nuclei to infer not only the presence of a closely separated black hole 
binary but also the mass ratio. When the black hole merger opens a hole in the inner disk, the broad band SED 
of the AGN or quasar may serve as a diagnostic. Such sources should be especially luminous in optical bands 
but intrinsically faint in X-rays (i.e., not merely obscured). We briefly note that viable candidates may have 
already been identified. 

Subject headings: black hole physics — galaxies: active 



1. INTRODUCTION 

The prevalence of supermassive black holes (SMBHs) at 
the centers of galaxies is well established (Richstone et al. 
1998). In the context of a hierarchical merging universe, 
merging galaxies can lead to the merging of the black holes 
(BHs Volonteri et al. 2003). Mergers of BHs are important as 
one pathway for the growth of BHs and have even been sug- 
gested as the principal cause of the scaling relations between 
BH mass and host galaxy properties (Jahnke & Maccio 201 1). 
Mergers are also strong gravitational wave emitters, and the 
asymmetric emission of gravitational waves, especially from 
spinning BHs, can lead to a recoil of the merged BH large 
enough to kick it out of the host galaxy potentially influencing 
the BH occupation fraction of galaxies (e.g., Volonteri et al. 
2008; Baker et al. 2008; Merritt et al. 2004). Even if two 
galaxies with BHs merge and the BHs sink to the center of 
the merged galaxy, however, it is not a given that the BHs will 
merge within a Hubble time (e.g., Milosavljevic & Merritt 
2003). 

The search for precursors to BH mergers has naturally fo- 
cused on active galactic nuclei (AGNs) that are either spatially 
resolvable at separations of ~ 1 kpc (e.g., Comerford et al. 
2012; Blechaetal. 2012) or have spectroscopically dis- 
tinct broad line regions at separations of ^0.1 pc (e.g., 
Gaskell 1996; Boroson & Lauer 2009; Smith et al. 2010; 
Eracleous et al. 201 1). 

One avenue for finding BH pairs at very small separations 
comes from an analogy with protoplanetary disks in which the 
presence of a planet may be inferred from its influence on the 
protoplanetary disk, particularly through gaps and holes in the 
disk carved out by the planet (Goldreich & Tremaine 1980; 
Takeuchi et al. 1996; Armitage & Rice 2005; Espaillat et al. 
2008). In this letter, we consider the observable consequences 
of an accretion disk gap caused by a secondary BH. Owing to 
the complexity of AGN and quasar spectra, we have devel- 
oped observational diagnostics in terms of powerlaw forms 
that can be fit to line-free regions in spectra, or potentially 



even to fluxes derived from sequences of broad filters. We 
have also given consideration to signatures that may be evi- 
dent in broad-band SEDs of quasars and AGN that span the 
optical and X-ray regimes. 

2. DESCRIPTION OF THE PHYSICAL SYSTEM 

The presence of a (proto-)planet in a proto-planetary disk 
in orbit around a central star can lead to transfer of angular 
momentum such that the planet moves inwards (Type I mi- 
gration; Goldreich & Tremaine 1980). Under certain circum- 
stances (section 3), the planet can clear material near its ra- 
dius, opening up a gap. Dynamical interactions between the 
planet and material at the gap edges tend to drive both the 
gap and the planet inward (type II migration; Takeuchi et al. 
1996). 

We are interested in the observational consequences of 
analogous systems for SMBHs. Consider a central BH of 
mass M with an accretion disk around it with a secondary BH 
of mass m, mass ratio q = m/M, in the plane of the disk at an 
orbital radius of a. If the secondary BH can open a gap in the 
accretion disk, the relative difference in emission would cre- 
ate an observational signature that could be used as evidence 
of close SMBH binary pairs at separations of ~ 10~ 3 pc. 

We assume that the accretion disk is as described by 
Shakura & Sunyaev (1973), i.e., an a-disk. The a-disk 
model is necessarily a mathematically convenient approxima- 
tion rather than a fully physical model, but it is still a suc- 
cessful model in many respects for a wide variety of astro- 
physical accretion systems. Our results are not bound to the 
Shakura & Sunyaev (1973) disk description, and we consider 
variations in section 5. What we require is that (i) azimuthally 
averaged disk midplane temperature is monotonic with radius 
in the inner ~ 10 3 R<j, where Rq = GMc~ 2 and (ii) h/r, the 
ratio of disk scale height to radius, is roughly constant. The 
recent considerations of azimuthal variations in disk temper- 
ature (Dexter & Agol 201 1; Dexter & Quataert 2012) still as- 
sume that the scale of the variability does not change with 
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radius so that for a given variability model, each azimuthally 
averaged annular emission is still a function only of r. The 
secondary hole is assumed to have already found its way into 
the accretion disk of the primary and settled into a coplanar 
orbit as has been shown to occur (Bogdanovic et al. 2007; 
Dotti et al. 2010). As we discuss below, the mass ratios of 
interest are q > 0.003. Below about q = 0.1-0.01, a poten- 
tial concern is that the mergers with halos containing small 
blck holes will lead to a complete dispersal of the secondary 
halo before the BH can settle to the center. Recent numeri- 
cal simulations, however show that it is possible for mergers 
with dwarf satellites can lead to secondary BHs of mass ratios 
q = 0.01-0.001 getting close enough to the primary that local 
processes can bring the BHs together (Bellovary et al. 2010). 
It is also possible that intermediate-mass BHs could form in 
situ in the accretion disk and would naturally be coplanar and 
have small mass ratios (McKernan et al. 2012). 

3. CONDITIONS FOR GAP OPENING 

There are two requirements that must be met in order for 
a gap in the disk to open: (1) the Hill sphere radius of the 
secondary must be larger than the disk scale height and (2) 
gap closing timescale from the viscous reaction of the disk is 
longer than the gap opening timescale. 

3.1. Hill Sphere Calculation 

The Hill sphere defines the region around an orbiting body 
in which it is gravitationally dominant and has radius: 
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The secondary BH will clear a gap in the accretion disk if 
Rh > h, which happens only when q>3(h/r) 3 . In the gas 
pressure dominated region where the binary orbits are most 
detectable, h/r = 0.004-0.008 for thin disks. We take a con- 
servative value of h/r = 0.01 as our fiducial value implies that 
q > 3 x 10~ 6 and even minor mergers are sufficient to carve a 
gap of width w ~ 2Rh- 

3.2. Gap closing 

Gaps will only exist if the gravitational torques open them 
faster than viscous diffusion can close them. The timescale 
for opening a gap 
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where Q. is the angular frequency of the secondary and m = 
rQ.c~ is the order of the Linblad resonance driving the gap 
creation (Takeuchi et al. 1996). The timescale to close the gap 
is 

„,2 



^close 



(3) 



where v is the kinematic viscosity. Requiring f c i ose > t open , 
yields the requirement that 
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where i2, c s , and h are all taken to be at the radius of the 
secondary. This requirement is equivalent to q > 10~ 5 for our 
fiducial disk parameters. 



3.3. Longevity of the gap 

In order for a gapped accretion disk to be of any observa- 
tional consequence, it needs to persist for a reasonably large 
fraction of an AGN's lifetime. If the angular momentum 
transferred from the secondary to the outer disk is greater than 
the angular momentum transferred from the inner disk to the 
secondary, then the secondary and the gap will both migrate 
inwards. The gap migration timescale when the local disk 
mass dominates over the secondary mass is (Takeuchi et al. 
1996) 
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For a 10 8 M e primary with a secondary at ~ IQ^Rq, this corre- 
sponds to ~ 3 x 10 5 yr. For higher mass secondaries, the gap 
migration time is lengthened in the satellite-dominated regime 
by a factor of a couple over this value so that the lifetime of 
the gap is ~ 10 6 yr. 

When the BHs are close to each other, gravitational ra- 
diation is strong enough to shrink the orbit on interesting 
timescales. The Peters (1964) orbit-averaged evolution of the 
semimajor axis of two masses in orbit around each other with 
eccentricity e is 
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where G is the gravitational constant and c is the speed of 
light. For the case we consider, e = 0. The width of the gap 
is Rh, and the timescale for the secondary BH to move out of 
the gap is 
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Taking a 
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and q = 10 2 , then T, 



2 x 10 7 (M/10 8 M. ) yr, much longer than other relevant 
timescales. 

So for a gap lifetime of ~ 10 6 yr and a typical AGN lifetime 
of ~ 10 7 yr (Martini 2004), we can see 10 _1 of AGNs with 
companions with the right mass ratio and with separations r < 
2 x 10 3 R G . 

4. EVOLUTION OF THE GAP 

Once a gap is established it will evolve along three possible 
routes when considering only viscous effects. At the begin- 
ning it is a gapped disk, effectively a standard disk with an 
annulus of width w missing. If the gap migration timescale is 
faster than or comparable to the timescale for accretion of the 
disk material inside of the gap, then the gap will have constant 
w/h with a decreasing radius. 

If the mass accretion rate for the inner disk is faster than 
the gap migration rate, then it will partially decouple from 
the outer disk. This will first lead to an effective widening of 
the gap from the inner edge, making it more detectable. The 
inner disk will be continuously depleted and replaced by any 
material that accretes across the gap. If the mass accretion rate 
across the gap is the same as if there were no gap, then there 
would be no difference. On the other hand, if the accretion 
rate is a small fraction of what it would have been otherwise, 
the inner disk would be replaced by an accretion disk of a 
smaller accretion rate. 
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The details of how much accretion will persist across a 
gap created by a satellite are the subject of current numeri- 
cal simulations. For the small mass ratios that are our pri- 
mary interest here, however, it is plausible that the accretion 
rate is only reduced to 10% of what it would be in the ab- 
sence of a satellite, i.e., m re d = 0.1m, where m = M /Mem and 
M Edd = LEddTT 1 ^ 2 «3x 10 N.V-/ M)M. yr" 1 for an as- 
sumed 77 =0.1. 

For higher mass ratios, it is plausible that very little matter 
will accrete past the gap, resulting in essentially no replace- 
ment of the inner gap. In such a case, instead of a gap, there 
would be a hole in central accretion disk. The observational 
consequences of an accretion disk with a central hole — and 
the possible electromagnetic signal once the BHs merge — 
have been discussed before (Milosavljevic & Phinney 2005; 
O'Neill et al. 2009), and we return to this in section 6. 

The time to accrete the entire inner disk is 

t a =M m M-\ (8) 

where M is the mass accretion rate and M m is the mass of the 
inner accretion disk. For our fiducial accretion disk, we take 
the BH to be accreting at a fraction m = 0.01. The largest 
orbit we consider is roughly at the boundary outside of which 
free-free opacity dominates over Thomson scattering so that 
the mass inside the gap is roughly 

M « 6 x 10- 11 I j£- J a- 4 /V 3 / 15 ,V7 (9) 
so that the accretion timescale is 

t a « 0.002a- 4 / 5 (jp) m 6/5 m 8 / 15 . (10) 

For our fiducial case, t a w 3 x 10 7 yr, so that the inner disk is 
depleting itself more slowly than the gap's inward micragion. 
The above is only valid for q <C 1 . For mass ratios near unity, 
the dynamical torques will clear a hole inside of the orbit on 
a short timescale. 

5. SPECTRAL ENERGY DISTRIBUTION OF A GAPPED ACCRETION 

DISK 

The spectral energy distribution (SED) of a thermally emit- 
ting disk with temperature a function only of radius is 

[■Rout 

F X (X) = B x [T{r')]g{r')2%r'dr', (11) 

JR m 

where B X (T) is the Planck function and g(f J ) is a function 
that describes the emissivity as a function of radius. We as- 
sume that g(r') = 1 (optically thick) where the disk material 
is present (i.e., R m < r 1 < r — w and r + w < r 1 < R OU [) and 
g(r) = elsewhere. 

The absence of a ring in a thermally emitting disk will show 
up in the SED as a broad, shallow dip. Without a gap, the 
SED will fall off as XF X ~ A -4 / 3 . With a gap, the SED wiU 
be appear roughly as a broken powerlaw with a slope that is 
steeper and shallower than A -4 / 3 blueward and redward of 
the maximum deviation from the standard SED. While it is 
obviously possible to model an observed SED with the proper 
Planck function treatment, we use a broken powerlaw model 
to demonstrate the precision required to identify a gapped 
disk. We generated SEDs (Fig. 1) for gapped disks using 
equation (11) for various values of w, r/Rc, M, and m. For 



each simulated SED, we fitted broken powerlaw in the vicin- 
ity of the break of the form 

(X~ x A<Ao 

where x > 4/3 and y < 4/3 with an arbitrary normalization 
implied. The change in slope, A = x — y, from one side of 
the break to the other at A = A) increases as the gap width 
increases. We find 

(13) 

where f(x) is a residual function that is smoothly varying over 
w/h and has absolute value \ f(x) — 1 1 < 0.35 for most of the 
parameter space of interest. The dependence on r, m, and M is 
easily understood since the peak wavelength of a blackbody 
scales as A ~ T and temperature and mass scale as T ~ 
TO V 4 r - 3 / 4 in the outer regions where the secondary spends 
most of the time and M ~ T 1 /4 . The dependence on w comes 
from the area weighting of the output SED. The change in 
slope is well approximated by 

A w 0.002(w//i) 5 / 4 (14) 

for w/h < 160. Above w/h > 160, the change in slope is 
larger than the above approximation because the emitting an- 
nulus just outside of the gap is further down the Wien expo- 
nential drop off. Since w/h is purely a function of q, it is pos- 
sible to infer the mass ratio of the binary from just the change 
in slope, independent of m, M, and r. 

Rafikov (2012) found that a gapped accretion disk would 
have an SED such that XF X « A 12 / 7 . If T(r) - r~ - 9 , then 
the peak of the SED scales as XF X « A~' 7 , approximately 
as A 12 / 7 . There is no qualitative difference in the approach to 
modeling SEDs of such gapped disks. We find that the change 
in powerlaw slope for such a temperature power-law is only 
decreased by about 10% compared to our fiducial tempera- 
ture profile. Thus one could use the overall slope of the con- 
tinuum powerlaw to identify gapped accretion disks and then 
use detailed measurements of A and A) to infer properties of 
the BHs. 

6. SED OF AN ACCRETION DISK WITH A HOLE 

For higher mass ratios and suitably small separations, al- 
most the entire region inside of the orbit will be dynamically 
unstable as can be seen from considering the size of the Roche 
lobes for q w 1 binaries. Thus all gas will be cleared out of 
this region, leaving an accretion disk with a hole of size r ~ 2a 
in the center. The effect of this on the SED is much more dra- 
matic than that of a gap. Most or all of the shortest wavelength 
emission will be cut out, leaving a red SED. 

The absence of an inner accretion disk owing to a BH 
merger event may therefore be identified through the broad- 
band SED. The outer accretion disk will still have a high ac- 
cretion rate, and should be luminous in optical bands. In- 
deed, the outer disk may be abnormally luminous in optical, 
or even in UV bands (see, e.g., Kocsis et al. 2012), depending 
on the mass accretion rate and the truncation radius. On sim- 
ple grounds, UV luminosity is not a good diagnostic: a holed 
disk around a BH of lower mass is observationally degenerate 
with a filled inner disk around a BH of higher mass. Rather, 
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Figure 1. All panels show simulated SEDs with gaps, with our fiducial parameters of M = 10 s M Q , h/r = 0.01, and assuming our standard temperature profile. 
The left and center panels assume r = IOOOSq and show the effect of increasing the gap width with values w/h = 0, 10, 60, 80, 100, 140, and 180. The central 
panel displays a zoomed in version of the same and shows the depth of of the dip in the SED increasing as well as the differences in slopes on either side of the 
local minimum increase as w/h increases. Empirically, we find that the change in slope scales approximately as (w//i) 5 / 4 , allowing one to infer the mass ratio of 
the binary from just the change in slopes. The right panel assumes w/h =180 and shows the effect of increasing the gap radius with values from r/Ro = 200 to 
2000 in steps of 200. For reference, also plotted is a gapless SED. Modeled as a broken powerlaw near the gap, we find that the break occurs at a wavelength that 
scales with gap radius as A ~ (r/Rg) 3 '' 4 . This comes from the temperature profile's dependence of T ~ r~ 3 / 4 and having the peak wavelength of a blackbody 
scale as A ~ T~' . The break wavelength also depends on M and m as shown in equation (13). 



diminished X-ray emission is a better diagnostic, since non- 
thermal X-ray emission is thought to tap into the deep grav- 
itational potential close to the BH. Residual accretion within 
the hole may persist at a low level, but even this gives rise to 
X-ray variability in Sgr A*. Thus, X-ray variability may be an 
additional hallmark. Thus, a source that is luminous in optical 
bands but intrinsically faint (not faint due to obscuration) and 
potentially variable in X-rays might signal a recent or pending 
BH merger event. 

It is possible that such sources have already been identi- 
fied. When the distance to a given AGN is known a priori, 
e.g., through spectroscopy, SDSS colors are able to distin- 
guish whether or not a source is intrinsically luminous in, 
e.g., red bands, or just apparently red through obscuration. 
In particular, g — i is very effective, or A(g — i), where one 
measures color relative to the mean quasar color in a given 
red-shift bin. Hall et al. (2006) employed this technique to 
3814 SDSS quasars with the aim of isolating obscured BHs. 
Follow-up observations with Chandra revealed that a sub- 
set of intrinsically red quasars are not X-ray-faint because 
of obscuration. In particular, seven sources show no evi- 
dence of even moderate intrinsic absorption in X-rays, and 
two are constrained to have internal column densities below 
2 x 10 21 cm" 2 (Hall et al. 2006). This suggests that such 
sources are rare (7/3814 = 1.8 x 10~ 3 ), consistent with ex- 
pectations based on the peculiar configuration required and 
the relevant timescales, and also consistent with expectations 
based on broader surveys (Gibson et al. 2008). Another inter- 
esting case for a quasar with a small inner hole in its accretion 
disk may be PHL 1811, which is optically very luminous but 
X-ray-faint (Leighly et al. 2004; Choi et al. 2005). 

7. DISCUSSION AND CONCLUSIONS 

We have described the basic observational appearance of a 
BH with a gapped or holed accretion disk, both of which are 
observable through their broad-band SEDs. We have focused 
our consideration on parameter space that results in observ- 
able signatures between 2000 Aand 2 \im. At shorter wave- 
lengths, absorption of ultraviolet photons will prevent unam- 
biguous determination of details of the AGN continuum. At 
longer wavelengths, reprocessed emission from dust grains 
will similarly contaminate thermal emission from an accre- 
tion disk. In this range, however, it is possible to measure the 
continuum emission to the precision needed to infer the pres- 



ence of a gap. Current and future large surveys such as SDSS, 
Pan-Starrs, and LSST may be able to exploit the observational 
signatures we have developed. 

This treatment has ignored some potential complications 
that will need to be considered in the future. Of chief con- 
cern is that we have treated the emission as only coming from 
the top of the disk. Emission coming from the walls of the 
gap will complicate the signal and, depending on the vertical 
temperature structure, either increase or decrease the observ- 
ability of the gap. Tidal features produced by the secondary 
as seen in simulations of protoplanetary disks will also need 
to be considered and may lead to additional corroborating ob- 
servational signals. The periods of the secondaries that we are 
observable are months to several years, allowing for a poten- 
tial detection of a modulating signal. 
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